Stripe phases - possible ground state of the high-T c superconductors 
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Based on the mean-field method applied either to the extended single-band Hubbard model or 
to the single-band Peierls-Hubbard Hamiltonian we study the stability of both site-centered and 
bond-centered charge domain walls. The difference in energy between these phases is found to be 
small. Therefore, moderate perturbations to the pure Hubbard model, such as next nearest hopping, 
lattice anisotropy, or coupling to the lattice, induce phase transitions, shown in the corresponding 
phase diagrams. In addition, we determine for stable phases charge and magnetization densities, 
double occupancy, kinetic and magnetic energies, and investigate the role of a finite electron-lattice 
coupling. We also review experimental signatures of stripes in the superconducting copper oxides. 

PACS numbers: 71.10.Fd, 71.27.+a, 74.25.-q, 74.72.-h 



INTRODUCTION 

Since the discovery of high-temperature superconduc- 
tivity by Bednorz and Miiller 0], the unusual physical 
properties of the copper oxides have stimulated theorists 
and have led to the appearance of many new ideas 0]. 
One of the especially appealing new pictures that has 
emerged is the instability towards a novel type of co- 
existing incommensurate (IC) charge and magnetic or- 
der, i.e., stripe phase. As a rare event in the theory 
of high temperature superconductivity, the theory pro- 
ceeded here the experiment and the existence of stripe 
phases was predicted on the basis of Hartree-Fock (HF) 
calculations in the two-band model for CuCb planes of 
layered La2- sSr^CuCU (LSCO) 0, before their exper- 
imental confirmation. This instability persists as well 
in the effective single-band Hubbard model 0, 000]. 
All these calculations yielded solutions with a phase sep- 
aration manifested in formation of nonmagnetic lines 
of holes, one-dimensional (ID) domain walls or stripes, 
which separate antiferromagnetic (AF) domains of oppo- 
site phases. Such states result from the competition be- 
tween the superexchange interaction, which stabilize the 
AF long-range order in the parent Mott insulator, and 
the kinetic energy of doped holes. Indeed, the magnetic 
energy is gained when electrons occupy the neighboring 
sites and their spins order as in the Neel state, whereas 
the kinetic energy is gained when the holes can move and 
the AF order is locally suppressed along a domain wall 
(DW). Thus, a stripe phase provides the best compromise 
between the superexchange promoting the AF order and 
the kinetic energy of doped holes. 

However, the debate on the microscopic origin of the 
stripe instability is far from closed. Two main scenarios, 
based on a Ginzburg-Landau free energy, for the driv- 
ingmechanism of the stripe phase have been discussed 
[abl- In the first one, stripes are charge-density waves 
with large periodicity arising from the Fermi surface (FS) 



instability with the transition being spin driven (j]. A 
general feature of such an instability is a gap / pseudogap 
which opens up precisely on the FS. Hence, the spacing 
between DWs is equal to 1/x, with x denoting doping 
level so as to maintain a gap/pseudogap on the FS. In 
this scenario spin and charge order occur at the same 
temperature or charge stripe order sets in only after spin 
order has developed. 

An alternative scenario comes from the Coulomb- 
frustrated phase separation suggesting that stripe for- 
mation is charge driven. Indeed, using the Ising model, 
it has been shown that the competition between long 
range Coulomb interactions and short range attraction 
between holes leads to formation of stripes [HJ- In this 
case Ginzburg-Landau considerations lead to an onset of 
charge order prior to spin order as the temperature is low- 
ered. However, the above analysis does not take into ac- 
count spin fluctuations which might be crucial for the na- 
ture of the phase transition by precluding the spins from 
ordering at the charge-order temperature [ijj. Moreover, 
the conjecture that long range Coulomb forces are re- 
quired to stabilize stripe phases has been challenged by 
the studies of the t-J model, in which the DW structures 
were obtained without such interactions (l^ . 

In order to investigate the influence of strong electron 
correlations due to large on-site Coulomb repulsion U at 
Cu ions, several methods have been employed to study 
the stripe phases which go beyond the HF approxima- 
tion, such as: Density Matrix Renormalization Group 
(DMR G) [H E3, Slave-Boson Approximation (SB A) 
[14Lll5llfij. variational local ansatz approximation [l7j . 
Exact Diagonalization (ED) of finite clusters ana- 
lytical approach based on variational trial wave function 
within the string picture [lflj . Dynamical Mean Field 
Theory (DMFT) HJ, Cluster Perturbation Theory 
(CPT) 0, and Quantum Monte Carlo (QMC) [HEl]. 
In spite of this huge effort, it remains unclear whether 
DWs are centered on rows of metal atoms, hereafter 
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named site-centered (SC) stripes, or if they are centered 
on rows of oxygen atoms bridging the two neighboring 
metal sites, the so-called bond-centered (BC) stripes, and 
even calculations performed on larger clusters did not 
yield a definite answer [2^]. Therefore, the purpose of this 
paper is to study the stability of both structures based 
on the mean-field method applied either to the extended 
single-band Hubbard model or the single-band Peierls- 
Hubbard Hamiltonian which includes the so-called static 
phonons [13] ■ For stable phases we determine charge 
and magnetization densities, double occupancy, kinetic 
and magnetic energies, and investigate the role of a finite 
electron-lattice coupling. 

EXPERIMENTAL SIGNATURES OF STRIPES 

Experimentally, stripe phases are most clearly detected 
in insulating compounds with a static stripe order, but 
there is growing evidence of fluctuating stripe correla- 
tions in metallic and superconducting materials. The 
most direct evidence for stripe phases in doped anti- 
ferromagnets has come from neutron scattering studies 
in which charge and spin modulations are identified by 
the appearance of some IC Bragg peaks, in addition to 
those which correspond to the crystal structure. How- 
ever, sometimes sufficiently large crystals are not avail- 
able for such experiments, and one has to resort to 
other methods capable of probing local order. These 
methods include nuclear magnetic resonance (NMR) , nu- 
clear quadruple resonance (NQR), muon spin rotation 
(/iSR), scanning tunneling microscopy (STM), and trans- 
mission electron microscopy (TEM). Furthermore, angle- 
resolved photoemission spectroscopy (ARPES), angle- 
integrated photoemission spectroscopy (AIPES), as well 
as x-ray photoemission (XPS) and ultraviolet photoemis- 
sion (UPS) spectroscopies all provide essential informa- 
tion about conspicuous changes in the electronic struc- 
ture when stripe structure sets in. Finally, a distinct 
imprint of the ID spin-charge modulation on transport 
properties should be detectable as the in-plane anisotropy 
of the resistivity and the Hall coefficient Rh ■ 

The abundance of the current evidence on various 
types of stripe order as well as the recent ARPES re- 
sults on the spectral weight of the cuprate superconduc- 
tors is contained in the review articles by Kivelson et al. 
[27^ . and by Damascelli et al. Historically, the first 

compelling evidence for both magnetic and charge or- 
der in the cuprates was accomplished in a neodymium 
codoped compound I^-z-yNdySr^CuC^ (Nd-LSCO). 
For y = 0.4 and x = 0.12, Tranquada et al. [H |s3 
found that the magnetic scattering is not characterized 
by the two-dimensional (2D) AF wave vector (1/2, 1/2), 
but by IC peaks at the wave vectors (1/2 ± e, 1/2) with 
e = 0.118. Moreover, inspired by the pioneering works 
demonstrating that the staggered magnetization under- 



goes a phase shift of it at the charge DWs 0, 0, [H0-0I- 
the authors found additional charge order peaks (±2e, 0), 
precisely at the expected position 2e = 0.236. Interest- 
ingly, this doping corresponds to a local minimum in the 
doping dependence of the superconducting temperature 
T c in Nd-LSCO suggesting that the static stripes 
are responsible for this anomalous depression of super- 
conductivity. However, it may well be that the apparent 
correlation is entirely accidental and therefore the role 
of stripes in superconductivity remains an open question 

Q 

Unfortunately, in early studies Tranquada et al. 
detected only magnetic IC peaks at higher doping lev- 
els x = 0.15 and x — 0.2. Nevertheless, systematic 
NQR studies of Nd-LSCO revealed the presence of robust 
charge stripe order throughout the entire superconduct- 
ing regime of doping 0.07 < x < 0.25 jsj. Also in a more 
recent study, both charge and spin superlattice peaks at 
x = 0.15 were found recently in the neutron diffraction 
experiments by Wakimoto et al. 

In fact, the reason why static stripes could be detected 
in this compound is a structural transition from the low 
temperature orthorhombic (LTO) to the low temperature 
tetragonal (LTT) phase, induced by the substitution for 
La ions by isovalent Nd ions. This, in turn, provides a 
pinning potential for dynamic stripes and stabilizes the 
charge order. Evidence of a similar pinning potential has 
also been found both in the /iSR and NQR studies of 
La 2 - a; _ y Eu !y Sr a; Cu04 (Eu-LSCO) with y ~ 0.2 |HE 
Moreover, the connection between the LTT phase and 
the appearance of charge and spin stripe order has been 
clearly demonstrated both in the neutron scattering and 
a;-ray diffraction studies on La2- a; -j,Baj,Sr a ;Cu04 (Ba- 
LSCO) with y = 1/8 [HE3|. Finally, static IC charge 
(2±2e, 0) and magnetic (l/2±e, 1/2) peaks have been de- 
tected within the LTT phase of La 2 -xBa a; Cu04 (LBCO) 
with x = 1/8 39] . The position of the peaks and the 
established incommensurability e = 0.118 are exactly the 
same as those obtained by Tranquada et al. for Nd- 
LSCO. Notably, the peaks that correspond to charge or- 
der appear always at somewhat higher temperature than 
the magnetic ones, indicating that the stripe order is 
driven by the charge instability. 

Let us now discuss the experimental evidence of slowly 
fluctuating stripes in La2_3;Sr K Cu04. The main differ- 
ence between the Ba and Sr codoped system is the fact 
that the latter undergoes a structural phase transition 
from the high-temperature tetragonal (HTT) phase to 
the LTO phase. As a consequence, in the superconduct- 
ing regime x > 0.06, the LSCO system exhibits purely dy- 
namic magnetic correlations which give rise to IC peaks 
at the wave vector (1/2 ± e, 1/2) specified in tetragonal 
lattice units 2-K/atetra- In their seminal inelastic neu- 
tron scattering studies, Yamada et al. established a 
remarkably simple relation e ~ x for 0.06 < x < 0.12, 
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FIG. 1: Summary of experimental data illustrating the dop- 
ing dependence of incommensurability e in the cuprates. Re- 
sults have been obtained by different groups : Nd-LSCO (Refs. 

EHMJpEl); LSCO (Refs. EE IS El III El LCO 
(ReElHl); Zn-LSCO (Refs. EllHl); YBCO (Refs^illll)- 

In LSCO, e has been denned as a distance from the IC peak 
position to the AF wave vector (1/2, 1/2) either in the or- 
thorhombic (x < 0.06) or tetragonal (x > 0.06) notation (see 
insets), whereas at x — 0.06, both definitions are used due to 
the coexistence of diagonal and parallel to the Cu-0 bonds 
spin modulations. 

followed by a lock- in effect at e ~ 1/8 for larger x. 

In contrast, in the insulating spin-glass regime of LSCO 
x < 0.06, quasielastic neutron scattering experiments 
with the main weight at zero frequency demonstrate 
that IC magnetic peaks are located at the wave vec- 
tors (1/2 ± e/V2, 1/2 ± e/V2) El El 0. This phe- 
nomenon has often been interpreted as the existence of 
static diagonal stripes, even though no signatures of a 
charge modulation were observed. Another possible ex- 
planation is the formation of a short ranged spiral or- 
der as its chirality also breaks the translational symme- 
try of the square lattice by a clockwise or anticlockwise 
twist E3 ■ Remarkably, even though the spin modulation 
changes from a diagonal to vertical/horizontal one, i.e., 
along Cu-0 bonds, at x around 0.06, e follows the doping 
x reasonably well over the entire range 0.03 < x < 0.12, 
as shown in Fig. ^ In fact, just for x = 0.06, both di- 
agonal (e = 0.053) and vertical/horizontal (e = 0.049) 
IC spin modulations have been found to coexist EH- In 
a stripe model this corresponds to a constant density of 
0.5 (0.7) holes per Cu atom in the DWs in the verti- 
cal/horizontal (diagonal) stripe phases, respectively, be- 
cause of the difference in Cu spacings in the two geome- 
tries, i.e., a or tho = V^atetra- In contrast, in the narrow 
region 0.02 < x < 0.024, IC magnetic peaks are located 
at the wave vector (1/2 ± e/2, 1/2 ± e/2) with e ~ x cor- 
responding to a constant charge of one hole/Cu ion along 



a diagonal DW 

EE EE E3i 

However, below x = 0.02, 
this does not hold anymore and the incommensurability 
gets locked with the value e ~ 0.014. 

Unfortunately, any concomitant charge ordering has 
not yet been detected in LSCO. Nevertheless, by com- 
paring the data based on the wipeout effect of 63 Cu 
NQR charge order parameter in LSCO with the ones 
obtained from charge stripe compounds as (Nd,Eu,Ba)- 
LSCO, Hunt et al. [48] concluded that a similar stripe 
instability exists in LSCO over the whole underdoped su- 
perconducting region 1/16 < x < 1/8. It is also worth 
mentioning that a very compiling evidence for its exis- 
tence has been established in the measurements of the 
in-plane resistivity and the dynamical infrared conduc- 
tivity anisotropy [iflllHrT]. 

Experimental detection of IC magnetic peaks in the 
LTO phase of LSCO suggests that the LTT structure is 
not essential for the appearance of stripes. This conjec- 
ture has been confirmed in experiments on the oxygen 
doped La2Cu04+5 (LCO) with the orthorhombic crystal 
structure (HJ. It is also supported by the evidence for 
static IC magnetic peaks in another orthorhombic com- 
pound La2-£ C Sr. r Cui_j,Znj,04 (Zn-LSCO) with y up to 
0.03, even though attempts to observe the charge order 
peaks were unsuccessful [H3,Ii3|- In fact, Zn substitution 
pins the stripe fluctuations similarly to the rare-earth el- 
ements. However, in contrast to the latter, it does not 
induce a structural transition to the LTT phase, but pro- 
vides randomly distributed pinning centers that promote 
meandering of stripes and correspondingly broadens IC 
peaks. 

An important question is whether charge stripes ap- 
pear solely in monolayered lanthanum compounds or if 
they are a generic feature of all the cuprates. The lat- 
ter conjecture seems to be supported by inelastic neu- 
tron scattering experiments on bilayered YBa2Cu306+a 
(YBCO) compounds that have identified the presence of 
IC spin fluctuations throughout its entire superconduct- 
ing regime In fact, as the doped charge is non- 
trivially distributed between the Cu02 planes and CuO 
chains, it is very difficult to determine the precise dop- 
ing level x in the Cu02 sheet of YBCO. Nevertheless, 
systematic studies by Dai et al. jH3| have shown that 
the incommensurability in YBCO increases initially with 
doping but it saturates faster than in LSCO, i.e., already 
at x ~ 0.1 with the value e ~ 0.1. Unfortunately, there 
is no any compelling explanation that would account for 
such a different behavior of e in both systems. Eventu- 
ally, charge order peaks have been observed in YBCO6.35 
but in spite of several attempts, no static charge order 
could be detected in YBCO6.5 and YBCO6.6 so far [55| . 

Furthermore, although some neutron scattering ex- 
periments have been performed on Bi2Sr2CaCu20s+5 
(BSCCO) sample, the sample has only produced weak 
evidence of the IC structure Ef|- In contrast, Fourier 
transform of the recent STM data has revealed some IC 
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peaks corresponding to a four-period modulation of the 
local density of states along the Cu-0 bond direction, 
which may imply the existence of stripes jH3|- Never- 
theless, definite answer pertinent to the appearance of 
stripes in all the cuprates remains still unsettled and fur- 
ther experiments are required to reach an unambiguous 
conclusion, even though the summary of the experimen- 
tal data illustrating the doping dependence of the incom- 
mensurability e in cuprates, depicted in Fig.^ includes 
an array of compounds. 

Tendency towards phase separation is also a starting 
point to understand the doping evolution of the electronic 
structure in LSCO and Nd-LSCO. For example, ARPES 
spectra measured at the X = (tt, 0) point in LSCO show 
that even though the data are solely characterized by 
a single high binding energy feature in the insulating 
regime, upon increasing doping one observes a systematic 
transfer of spectral weight from the high- to the low bind- 
ing energy part [H3|. Consequently, a well-defined quasi- 
particle (QP) peak develops near the optimal doping. In 
contrast, the intensity near the S = (tt/2,tt/2) point re- 
mains suppressed for the entire underdoped regime so 
that a QP peak is observed only for x > 0.15. 

Another peculiar feature of the ARPES band disper- 
sion is extensively discussed in the literature saddle point 
at the X point, the so-called flat band As hole dop- 
ing increases, the flat band moves monotonically upwards 
and crosses the Fermi level Ep at x ~ 0.2. This is re- 
flected in the enhancement of the DOS at the chemical 
potential N(fi) observed by AIPES 0. 

The experimental distribution of the photoemission 
spectral weight near the X and S points in doped LSCO 
has been nicely reproduced using the DMFT approach for 
vertical SC stripes obtained within the Hubbard model 
[icj . As a consequence of the stripe order, the obtained 
spectra along the r — X — M path were not equivalent 
to those along the T — Y — M one, with T — (0, 0) and 
Y = (0,7r). Moreover, as in the experiment, the spectral 
weight along the T — X direction was suppressed close 
to the r point and simultaneously enhanced at the X 
point. Furthermore, in the framework of stripes, the flat 
QP band near the X point with a large intensity at the 
maximum below the chemical potential fi follows from 
a superposition of the dispersionless ID metallic band 
along the x direction, formed by holes propagating along 
the vertical domain walls, and an insulating band that 
stems from the AF domains. In contrast, an AF band at 
the Y point is characterized by a high binding energy well 
below /k and consequently the spectral weight at oj = fi 
almost vanishes. Finally, a distinct gap for charge exci- 
tations should open at [x near the S point. This gap fol- 
lows indeed from the stripe structure — while the system 
may be metallic along the stripes, i.e., in the antinodal 
directions T — X or T — Y, the low-energy excitations 
should be noticeably suppressed along the nodal direc- 
tion F — S crossing all the stripes. This conjecture is 



also supported either by the ED studies l.S| or by the 
analytical approach based on variational trial wave func- 
tion within the string picture 0], both applied to the 
t-t'-t"-J model, or by the CPT for the t-J model [13]. 

In fact, the low-energy spectral weight of Nd-LSCO 
at x = 0.12, a model compound for which the evidence 
of spin and charge stripe order is the strongest, is also 
mostly concentrated in flat regions along the V — X and 
r — Y directions, while there is only little spectral weight 
along the T—S direction [filj . On the other hand, ARPES 
spectra of both LSCO and Nd-LSCO at x — 0.15 have 
revealed not only the presence of flat bands around the 
X and Y points, but also the existence of appreciable 
spectral weight at Ep in the nodal region jjj^. While the 
observation of flat seg ments might be directly ascribed to 
ID domain walls detection of nodal spectral weight 
poses a formidable task to develop a theory that would 
describe the electronic structure resembling the FS of a 
fully 2D system because, as it was already stressed out, 
the nodal spectral weight is expected to be suppressed 
in a static SC stripe picture [H, E3, El- Indeed, 
the experimentally established FS looks rather like the 
one arising from disorder or from dynamically fluctuating 
stripes |fi3j . 

Alternatively, guided by the CPT results showing that 
while the SC stripes yield little spectral weight near the 
nodal region, the BC ones reproduce quite well the nodal 
segments [22|, Zhou et al. 62] have conjectured that 
the experimental FS may result from the coexistence of 
the SC and BC stripes. Within this framework, upon in- 
creasing doping the BC stripes are formed at the expense 
of the SC ones. This scenario is particularly interesting 
because it has been shown that the BC stripe, in contrast 
to its SC counterpart, enhances superconducting pairing 
correlations The relevance of a bond order at the 

doping level x = 0.15 is supported by recent studies of 
the ARPES spectra in a system with the BC stripes [f?5l |. 
These studies have yielded pronounced spectral weight 
both in the nodal and antinodal directions, reproduc- 
ing quite well the experimental results in Nd-LSCO and 
LSCO |^|. Furthermore, the stripe scenario would also 
explain the origin of the already discussed two compo- 
nents seen in the ARPES spectra at the X point near 
x = 0.05 [H3|. Indeed, the response from the AF insulat- 
ing regions would be pushed to the high binding energies 
due to the Mott gap, whereas the charge stripes would 
be responsible for the other component near Ep. 

Existence of DWs should also give rise to the appear- 
ance of new states inside the charge-transfer gap that 
would suppress the shift of the chemical potential fi in 
the underdoped regime x < 1/8 where e increases lin- 
early. Such pinning of u in LSCO was indeed deduced 
from XPS experiments In contrast, in the overdoped 
region with a lock-in effect of e, the number of stripes per 
unit cell saturates, doped holes penetrate into the AF do- 
mains, and consequently /i would move fast with doping 
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in agreement with the experimental data. The picture of 
broadened stripes and holes spreading out all over the AF 
domains above x = 1/8 is also indicated by the doping 
dependence of the resistivity and the Hall coefficient Rh 
in Nd-LSCO. Namely, a rapid decrease in the magnitude 
of Rh for doping level x < 1/8 at low temperature pro- 
vides evidence for the ID charge transport, whereas for 
x > 1/8, relatively large Rh suggests a crossover from 
the ID to 2D charge transport [67). Altogether, it ap- 
pears that the metallic stripe picture does capture the 
essence of the low-lying physics for Nd-LSCO and LSCO 
systems. 

Conversely, it is important to note that so far no ev- 
idence of IC peaks has been detected in any electron- 
doped cuprates superconductors. Instead, the neutron 
scattering experiments have established only commen- 
surate spin fluctuations as in Nd2- £! ;Ce a ;Cu04 (NCCO), 
both in the superconducting and in normal state [68l |. 
Moreover, observation of such peaks is consistent with 
the XPS measurements in NCCO showing that the chem- 
ical potential increases monotonously with electron dop- 
ing 



NUMERICAL RESULTS 

In this Section we attempt a systematic investigation 
of the properties and relative stability of filled vertical 
and diagonal stripes. We shall see that in spite of the 
difficulty to stabilize the ground state with half-filled 
stripes (one hole per every two atoms in a DW), the 
mean-field framework is useful as providing a generic mi- 
croscopic description of filled inhomogeneous reference 
structures with the filling of one doped hole per stripe 
unit cell. Their special stability rests on a gap that opens 
in the symmetry broken state between the highest occu- 
pied state of the lower Hubbard band and the bottom 
of the so-called mid-gap bands, i.e., some additional un- 
occupied bands lying within the Mott-Hubbard gap that 
are formed due to holes propagating along DWs (26j . 

Here, we extend early HF studies of the filled DWs 
[j. LJ. lJ LjI and determine a phase diagram of the Hub- 
bard model with an anisotropic nearest-neighbor hopping 
t by varying the on-site Coulomb repulsion U and inves- 
tigating locally stable structures for representative hole 
doping levels x — 1/8 and x — 1/6. We also report 
the changes in stability of the stripe structures in the 
extended Hubbard model due to the next-neighbor hop- 
ping t' and to the nearest-neighbor Coulomb interaction 
V. Finally, in order to gain a comprehensive understand- 
ing of the competition between different types of stripes 
in a realistic model, we include lattice degrees of freedom 
induced by a static Peierls electron-lattice coupling. 



Extended single-band Hubbard model 

The starting point for the analysis of stripe structures 
is the extended single-band Hubbard model, which is 
widely accepted as the generic model for a microscopic 
description of the cuprate superconductors |xU, 
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where the operator c\ a (Cj a ) creates (annihilates) an 
electron with spin a on lattice site i (j), and rn — 
cLc iT + cJjCji gives the electron density. The hopping 



tij is t on the bonds connecting nearest neighbors sites 
and t' for second-neighbor sites, while the on-site 
and nearest-neighbor Coulomb interactions are, respec- 
tively, U and V. 

The model can be solved self-consistently in real space 
within the HF, where the interactions are decoupled into 
products of one-particle terms becoming effective mean 
fields that act on each electron with the same strength. 
This approximation basically involves solving an eigen- 
value problem. The obtained wavefunctions form a new 
potential and hence the Hamiltonian for a new eigenvalue 
problem. Typically, the new potential is chosen as some 
linear combination of the current and preceding poten- 
tial. The iterations are continued until the input and 
output charge density and energy do not change within 
some prescribed accuracy. The most significant drawback 
of this method is that it neglects correlations. Electron 
correlation changes the system properties and manifests 
itself in the decrease of the ground state energy. The dif- 
ference between the energy of the exact ground state and 
the energy obtained within the HF is thus called the cor- 
relation energy. It arises from the fact that an electron's 
movement is correlated with the electrons around it, and 
accounting for this effect lowers further the energy, be- 
yond the independent electron approximation. 

We do not consider noncollinear spin configurations, 
and use the most straightforward version of the HF with 
a product of two separate Slater determinants for up and 
down spins, whence, 



(2) 



A similar decoupling is performed for the nearest- 
neighbor Coulomb interaction. Calculations were per- 
formed on 12 x 12 (16 x 16) clusters for x = 1/6 (x = 1/8) 
with periodic boundary conditions, and we obtain stable 
stripe structures with AF domains of width five atoms 
for x = 1/6 and seven atoms for x = 1/8. Typical solu- 
tions at x = 1/8 are shown in Fig. with the local hole 
density, 



(n hi ) = 1 - (n iT + riii), 



(3) 



scaled by the diameter of the black circles and the length 
of the arrows being proportional to the amplitude of local 



>'-7->'-'- 



)t r ,'-^'-,^-7'-,'--^- 

^'777^ 



/////' 



^->7->7'^ 



\yLy'-y'-y<-,--,'- 
^7->7^- 

y-y'-y'-y'-,-^-- 



7^' 



///// 



>7? 



y'-y'-y'-y'-y'T- 



7-y'V-y'-7- 
///' 



> y 7-/77} 

7/7/7/7 



■ '-y--y'-y'-y'-y--y'- 



,<-,'-,'-,^-y'---^' 



7y'-y'V-y'-y'-y'—7 



JfryYi' 
±>hY?Y 



//////■' <•//' 

T-^-y'-yh'-y'-y' 

7y'-y'-y'-y--y'-y-- 



77 77 7 > I >» 




/ 
/ 

y 
/ 
y 
/ 
/ 
/ 
/ 
/ 
y 
/ 
y 
/ 
y 
> 

/ 
/ 

y 

/ 
/ 

y 
/ 

y 


h 




t-, 
f 




^y 


'-7 


^y 


Ly 


Ly 


h 

Ly 


* 







Ly 


Ly' 


Ly 








<-y 


<-, 






Ly 










Ly 




Ly 








-i 




Ly 


<h 

Ly 


Ly 


<h 


f 




y 


H 


% 
Ly' 










'-? 


'-y 


Ly 




Ly 


fy 






✓ 


Ly 




Ly 








^y 




Ly 


h 


Ly 


<h 

Ly 






y 


*-> 


V 

Ly' 






/ 




'-> 


'-y 




Ly 


Ly 

Ly 


Ly 


f 






Ly 




Ly 












Ly 


Ly 




Ly 


f-, 




y 


Ly 


Ly' 
Ly' 


Ly 


Ly 








'-y 


<K 


Ly 


*1 

Ly 




* 






Ly 




Ly 

Ly 








~y 


h 


Ly 


Ly 

Ly 




'-y 


f 




y 


^y 


>:< 




Ly 

'-, 






<-y 


'-y 


Ly 






'-y 


'—, 






Ly 




h 

Ly 








-> 


<-, 


Ly 


Ly 


Ly 


Ly 


/ 




y 


^y 


>' 




«7 

<-, 






'-> 


'-y 


Ly 




Ly 


<Ly 


t—, 




/ 


Ly 






Ly 










Ly 


<h 


Ly 


Ly 


t-, 




y 


fy 


t' 


Ly 


^ 








H 
'-, 


Ly 


Ly 


Ly 

Ly 


Ly 


t 




/ 


Ly 




Ly 








-y 






Ly 






I 




y 


Ly 
* 


% 


s / 
/ j . 


' r * / / t * 

* y / / t S J . 




/ 






t 


/ 




/ 

h 




* 








/ 




A' 






' / 


'-> 






Ly 


/ , 








Ly 




h 

Ly 






' f 


^7 




^y 


'-y 








/ 


Ly 


Ly 






Ly' 
/ / 


'/ 


h, 




h 




< 






/ 


Ly 


Ly 


h 


h 






Ly 






'-, 

'-, 


' , 






/ 

Ly 


Ly 


Ly 


•Ly 

Ly 


■-y 
'-y 


> 








h 




<-, 










Ly 


Ly 




' 






Ly' 




Ly 


<-, 




-y 


-> 


Ly 


Ly 


Ly 


Ly 


^y 

f 






* / 




*7 


/ f 






-y 




Ly 


Ly 


Ly 


y^ 








Ly 


^ 


Ly 

<-. 






/ 






*J 


Ly 
Ly 


*y 
/ 








Ly 










/ 




^y 


'-y 


Ly 

Ly 








' > 






Ly 


Ly' 

Ly' 


y 
/ 

/ 




/ 




h 




'-y 








' / 






Ly 


h 


'-, 


/ 










< 






' / 








Ly 

Ly 


Ly 




'/ 






h 










Ly 


Ly 






/ , 




^y' 


Ly 


h 










v t 


Ly 


ty 














y 


Ly 
/ 








f 


Ly 


h 


Ly 


H 


* 






/ 


Ly' 


Ix 




PIG. 2: Vertical site-centered (VSC) and diagonal site- 
centered (DSC) stripe phases as found for U/t = 5 at hole 
doping x = 1/8. The length of arrows is proportional to the 
magnetization (Sf) and the hole density (rihi) is scaled by the 
diameter of black circles. 



FIG. 3: Vertical bond-centered (VBC) and diagonal bond- 
centered (DBC) stripe phases as found for U/t = 5 at hole 
doping x = 1/8. The meaning of the arrows and black circles 
as in Fig. El 



magnetization density, 

<S?> = ±1(^-^)1. (4) 

These structures possess nonmagnetic DWs with en- 
hanced hole density which separate AF domains having 
hole density almost unchanged with respect to the un- 
doped case. Note that the AF sites on each side of the 
DWs have a phase shift of n. 

In order to appreciate better the microscopic reasons 
of such arrangement let us consider a small cluster con- 
sisting of three atoms filled by two electrons and one 
hole (with respect to half-filling with the electron density 
n = 1 per site). For simplicity we assume that the elec- 
trons are confined to the considered cluster owing to large 
Coulomb interaction U ^> t, and we do not take into ac- 
count any interactions with the AF background. There 
are two possible candidates for the ground state. The 
first one corresponds to a hole added to three atoms of a 
single AF domain in which, if we suppose that a |-spin 
electron is replaced by a hole, the two remaining f-spin 
electrons can be found in one of three allowed configura- 
tions: {f, 0, t}, {Ti T? 0}, and {0, |, f } (the other configu- 
rations are excluded by the Pauli principle). Hence, this 
polaronic state gives the total energy, 

E P = -V2t, (5) 

and the Coulomb interaction U does not contribute. 



A different situation is obtained when a hole occupies 
instead a DW separating two AF domains. Derealiza- 
tion leads then to similar three configurations to those 
obtained above: {T,0, J.}, {1\!,0}, and {0,t,|}, but in 
addition, three configurations with one doubly occupied 
site {Ti, 0, 0}, {0, tl, 0}, and {0, 0, T!}, can be reached as 
excited states which cost Coulomb energy U. Moreover, 
three other configurations with interchanged |- and |- 
spins are then also accessible via the decay of double oc- 
cupancies: {J.,0, T}, {|,T)0}, and {0, |,T}- In the regime 
of large U, the total energy in the ground state can be 
found in a perturbative way, and as a result one obtains, 



E s = -V2t-—. (6) 

Therefore, the Hilbert space for the latter solitonic solu- 
tion is larger and one finds that this solution is always 
more stable than the polaronic one j2^|. The argument 
applies also to 2D systems, where the DWs are more sta- 
ble than the lines of polarons in an AF background. 

We compare the stability of such nonmagnetic SC do- 
main walls with the BC stripe phases in which DWs are 
formed by pairs of magnetic atoms, as obtained by White 
and Scalapino 

EH (cf- Fig. EJ. In the three-band model, 
SC (BC) stripes correspond to DWs centered at metal 
(oxygen) sites, respectively [HI IzlIzS Izil • 
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Effect of hopping anisotropy 

We begin by setting t' = and V = with the goal 
of elucidating the effects of hopping anisotropy on the 
stripes. This is motivated by the fact that the first detec- 
tion of static stripes in both charge and spin sectors was 
accomplished in Nd-LSCO indicating that rare-earth 
elements doping is in some way helpful for pinning the 
stripe structure. Indeed, it produces a structural transi- 
tion in the system from the LTO to LTT phase (zjj. Both 
phases involve a distortion of the CuC>2 plane by rotation 
of the CuC>6 octahedra. In the LTO phase the tilt axis 
runs diagonally within the copper plane, such that all the 
oxygen atoms are displaced out of the plane. Conversely, 
in the LTT phase this rotation takes place around an axis 
oriented along the planar Cu— O bonds. Therefore, oxy- 
gen atoms on the tilt axis remain in the plane, while the 
ones in the perpendicular direction are displaced out of 
the plane. This provides a microscopic origin for in-plane 
anisotropies — the Cu— Cu hopping amplitude t depends 
on the Cu— O bond and it is isotropic in the LTO phase 
and anisotropic in the LTT one. For a physical tilt angle 
of order 5°, the relative anisotropy taking t y < t x , 



ty\ 



(7) 



is weak and amounts to et — 1.5% [zHIzS- The direc- 
tion with a larger hopping amplitude coincides with the 
direction of a stronger superexchange coupling J. 

The possible relationship between this anisotropy and 
the onset of stripe phases has been intensively studied 
within anisotropic Hubbard (t x ^ t y ) or t-J (t x ^ t y , 
Jx 7^ Jy) models by means of various techniques: un- 
restricted HF approach [tfj, DMRG 0, and QMC 
method j^. The in-plane anisotropies might also be 
represented theoretically by on-site potentials as in the 
QMC study by Riera |2J|. All these studies have shown 
a pronounced tendency to forming stripe phases, which 
manifests itself by the reduction of their energy |7jj,|7j|, 
accompanied by the appearance of IC peaks in the spin 
and charge structure factor 0|- K appears that 
a finite anisotropy of the next-nearest hopping term t' 
might play a role in stabilizing diagonal incommensu- 
rate peaks observed in the spinglass phase of LSCO 
(0.02 < x sj 0.06) mmSlialllQ. Indeed > al " 
though the LTO phase is usually considered as isotropic, 
which is the case for nearest-neighbor hopping and inter- 
action, a different length of the orthorhombic axes implies 
the need for an anisotropic t' parameter. Exact diagonal- 
ization studies incorporating such anisotropy have shown 
that it strongly strengthens hole correlations along one 
direction and suppresses them along the other, resulting 
in a ID pattern of holes [7s| . 

It turns out, however, that the variation of the hopping 
anisotropy e t Q has only a little visible effect on the local 
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FIG. 4: Local hole n h (l x ) (top) and magnetization S n Qx) 
(second row) density; kinetic energy Ef(l x ) (third row) 
and Et{l x ) (bottom) projected on the bonds in the ^(in- 
directions, respectively, of the VSC (left) and DSC (right) 
stripe phases shown in Fig. [5] (open circles) as well as of the 
ones obtained in the anisotropic model with t x /t y = 1.22 
(filled circles) . For clarity, the latter are shifted by one lattice 
constant from the origin of the coordinate system. 



hole density, 



n h (l x ) = 1 - (ra(i s ,o),T + ra^.ou): 



(8) 



shown in Fig. 01 as a function of the a;-direction coordi- 
nate l x for a given y-direction coordinate l v = 0, even at 
the unrealistically large anisotropy level e t = 22%, cor- 
responding to t x /t — 1.1 and ty/t = 0.9. Similarly, the 
anisotropy does not modify the modulated magnetization 
density, 



S w (l x ) = (-1)^1(^,0),? - n (i x ,o),i)> 



(9) 



with a site dependent factor (— l) ix compensating mod- 
ulation of the staggered magnetization density within a 
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FIG. 5: The same as in Fig.Q]but for the BC stripe phases. 

single AF domain. 

In contrast, the strong effect of finite anisotropy e* J2J 
is clearly demonstrated by variation of the expectation 
values of the bond hopping terms along the x- and y- 
directions, 

E?(l x ) = -t x Y,( c li^),. c (i*+i,o),* + h.c), (10) 
= -h Y,( c Lou c (LA),« + M- (11) 





Ef/t 


E y Jt 


Eu/t 


Etat/t 


VB(S)C 


-0.6753 


-0.6147 


0.4900 


-0.8000 


DBC 


-0.6375 


-0.6375 


0.4726 


-0.8024 


DSC 


-0.6368 


-0.6368 


0.4696 


-0.8040 



TABLE I: Site- normalized ground-state energy E to t, kinetic 
energy (Ef,E%), and potential energy Ejj in the isotropic 
Hubbard model with U/t = 5 and i = 1/8 as obtained for 
different stripe phases: vertical site-centered (VSC), diagonal 
site-centered (DSC), vertical bond-centered (VBC) and diag- 
onal bond-centered (DBC). In the HF, both types of vertical 
stripes are degenerate. 



the anisotropic model. Analogous conclusion based on 
Fig. might be drawn concerning the orientation of the 
VBC stripes. 

Regarding diagonal stripes, although a finite 
anisotropy in hopping is also reflected in the ki- 
netic energy anisotropy, a system with either the DSC 
or DBC stripe pattern becomes topologically frustrated 
and consequently may gain less kinetic energy compared 
to a system with vertical stripes, taking a full advantage 
of the hopping anisotropy ( cf. Tables III and ITTf) . 

The effect of an increasing anisotropy illustrates the 
phase diagram shown in Fig. U3 determined by varying U 
and the ratio t x /t y of the nearest-neighbor hoppings in 
the x- and y-directions, while maintaining constant t = 
\{t x +t y ). We observe the generic crossover from vertical 
to diagonal stripes with increasing Coulomb interaction 
reported in early HF studies 0, [j E3- 01- The transition 
from the VSC to DSC stripes appears in the isotropic 
case at U/t ~ 4.1 for x — 1/8, and at a higher value 
U/t~ 4.6 for x = 1/6 [cf. Fig. Eta)]. The corresponding 
phase boundary between the VBC and DBC stripes is 
shifted towards stronger Coulomb interaction and occurs 
at U/t ~ 4.4 (5.0) for x = 1/8 (x = 1/6), respectively 
[cf. Fig.Etb)]. 

The results shown in Fig. have a simple physical in- 
terpretation. Stripe phases occur as a compromise be- 
tween, on the one hand, the AF interactions between 
magnetic ions and the local Coulomb interactions which 
favor charge localization, and the kinetic energy of doped 
holes which favors charge derealization on the other 
hand. The kinetic energies in Table G] show further that 



These features are seen in Fig.Q] For the VSC stripes one 
finds a large anisotropy in the values of the kinetic ener- 
gies II10II and Ijllll . which becomes especially pronounced 
beside the stripes, and is strongly reinforced by the hop- 
ping anisotropy. Therefore, taking into account that the 
hopping between two different charge densities is favored 
over motion between equal densities, one should expect 
that transverse charge fluctuations will always tune the 
direction of DWs along the weaker hopping direction in 





Ef/t 


E//t 


Eu/t 


Etat/t 


DBC 


-0.8143 


-0.4807 


0.4815 


-0.8135 


DSC 


-0.8098 


-0.4836 


0.4793 


-0.8141 


VB(S)C 


-0.8304 


-0.4776 


0.4938 


-0.8142 



TABLE II: The same as in Table [H but with the hopping 
anisotropy e t = 22%. 
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FIG. 6: Phase diagrams for stable: (a) site-centered (SC), 
and (b) bond-centered (BC) stripe structures obtained in the 
anisotropic Hubbard model on a 16x16 cluster for doping 
x = 1/8 (solid lines) and on a 12x12 cluster for x = 1/6 
(dashed lines) Parameters: t' = 0, V = 0. 



the vertical stripes are more favorable for charge dynam- 
ics. This result, which is not immediately obvious, has 
however a straightforward origin. Namely, the HF always 
leads to a large spin polarization since it is the only way 
to minimize the on-site Coulomb repulsion. Indeed, re- 
moval of a !-spin electron at site i leads to relaxation of 
the f-spin electron energy level at this site. As a conse- 
quence, an alternating on-site level shift develops yielding 
an energetical motivation for the symmetry breaking and 
forming the AF order. 

However, the renormalization of the double-occupancy 
energy involves a strong reduction of the kinetic energy 
in the |-spin channel between site i and its neighboring 
sites, as an electron incoming into this site encounters 
a high energy potential U(n^). Therefore, in the HF 
approximation we shall be able to identify dynamically 
favorable stripe patterns only by comparing appropriate 
local magnetization densities. For example, charge fluc- 
tuations occur more readily in the VSC stripe geometry 
presumably due to their greater overall width indicat- 
ing weaker correlation effects (c/. Fig. 31 ■ This explains 
their stability at small U where the consequent cost in 
potential energy Ejj becomes insignificant. By contrast, 
the DSC stripes are narrower having larger hole density 
along nonmagnetic DWs. Moreover, magnetization den- 
sity of their nearest neighbor sites is markedly enhanced 
as compared to the corresponding VSC stripe magneti- 
zation, as shown in Fig. and in Table ITTT1 The former 
also illustrates that the bonds connecting DWs with their 
nearest neighboring sites perpendicularly to the walls, 
have the main contribution to the kinetic energy gain, in 
fact suppressed here by larger spin polarization. Taken 
together, the above features are reflected in a more lo- 
calized character of the DSC stripes, with a lower net 
double occupancy and hence a more favorable on-site en- 
ergy Ejj (c/. Table||J|. This clarifies the mechanism of the 
transition from the VSC to DSC stripes with increasing 
U. 

Turning now to the analogous crossover between the 



i 


1 2 3 4 5 


w 


0.364 0.234 0.067 0.014 0.006 


VSC 


(0.378) (0.234) (0.060) (0.013) (0.006) 


(Si) 


0.000 0.222 0.348 0.381 0.384 




(0.000) (0.234) (0.357) (0.382) (0.384) 


(n hi ) 


0.388 0.193 0.070 0.032 0.020 


DSC 


(0.405) (0.195) (0.066) (0.028) (0.017) 


(Si) 


0.000 0.262 0.352 0.373 0.380 




(0.000) (0.272) (0.360) (0.377) (0.382) 



TABLE III: Local hole (n hi ) and magnetization (Si) density 
of the site-centered stripes shown in Fig. [2] all labeled by 
decreasing hole density in the a;-direction. In parenthesis the 
values for the extended hopping model with t' ft — —0.15 are 
given. 

BC stripes, we shall again compare local hole and magne- 
tization densities on and around their DWs. In contrast 
to the SC case, a VBC stripe phase possesses larger hole 
density along DWs, as illustrated in Fig. Eland Table IWl 
suggesting that it is more localized than the DBC one. 
Nevertheless, a better renormalization of the double oc- 
cupancy energy Eu by the latter (c/. Table QJ follows 
from a stronger spin polarization not only of the DW 
atoms but also their nearest neighbors (cf. Fig. 03 and 
Table IWll . This enhancement is directly responsible for 
a substantial reduction of the kinetic energy along bonds 
joining these atoms. Correspondingly, it accounts for a 
crossover from the DBC to VBC stripes in the small U 
regime when the larger kinetic energy gain becomes cru- 
cial. 

We would like to emphasize that the above transition 
between different types of stripe phases is not an artefact 
of the HF and occurs also between filled stripes obtained 
within more realistic approaches including local electron 
correlations. Indeed, slave-boson studies of the Hubbard 
model at the doping x = 1/9 have established that the 
transition from the filled VSC to DSC stripe phase ap- 
pears at the value U/t ~ 5.7, being much higher than 
that predicted by the HF, which yields U/t ~ 3.8 Q. 



i 


12 3 4 


VBC (n hi ) 


0.326 0.136 0.030 0.007 


(S!) 


0.118 0.301 0.371 0.384 


("hi) 


0.314 0.115 0.047 0.023 


DBC 


(0.323) (0.110) (0.046) (0.021) 


(si) 


0.145 0.322 0.365 0.378 




(0.155) (0.333) (0.368) (0.380) 



TABLE IV: The same as in Table EH but for the bond- 
centered stripes. VBC stripe is unstable in the extended hop- 
ping model with t'/t = —0.15. 
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In this method, enhanced stability of the VSC stripes fol- 
lows from an additional variational parameter per each 
site di, reducing the on-site energy without a strong sup- 
pression of the kinetic energy. Remarkably, the total en- 
ergy difference between the vertical SC and BC stripes 
at both doping levels is comparable to the accuracy of 
the present calculation. Such degeneracy was also re- 
ported in the HF studies of the charge-transfer model 
|71| . However, when electron correlations are explicitly 
included the BC stripes are more stable at and above 
at = 1/8 doping [Ulza]. 



Effect of the next-neighbor hopping t' 

We now turn to the effect of a next-neighbor hopping if 
on the relative stability of the stripes. There are numer- 
ous experimental and theoretical results which support 
the presence of finite t' in the cuprates. For example, re- 
cent slave-boson studies have revealed that the phenom- 
ena of the half-filled vertical stripes in LSCO requires a 
finite next-neighbor hopping t'/t ~ —0.2 [lfj . 

Let us pause now for a moment to clarify the influence 
of t' on the DOS as well as on the FS using the electronic 
band which follows from a simple tight-binding model, 

E(k) — — 2i(cos k x + cos k v ) — 4t' cos k x cos k y . (12) 

By the reduction from the Cu02 multiband model to an 
effective single-band model it has been found that t > 
and t' < for hole doped system, and t < and if > in 
electron doped system [ziJ]. Although an accidental can- 
cellation of the various contributions results in almost 
perfect electron-hole symmetry of the nearest-neighbor 
hopping t, the next-neighbor hopping t' asymmetry ap- 
pears owing to the fact that the dominant contribution 
to the latter comes from a direct 0-0 hopping t pp in the 
case of a hole hopping. On the contrary, an electron hop- 
ping follows from a third order Cu^O^O^Cu process, 
being therefore dominated by the Cu-0 hopping element 

ipd' 

In the noninteracting limit the role of if is to shift the 
van Hove singularity away from the middle of the band, 
either to higher or to lower energy depending on its sign 
[scj . Fig. shows the tight-binding DOS, centered at 
lo = with the condition J N(oj)u)dijj = 0, and the oc- 
cupied states at the doping x = 1/4. In the hole-doped 
case, with the vacuum as the zero electron state, the van 
Hove singularity lies in the lower part of the band. Con- 
versely, in the case of electron doping, with the vacuum 
as the zero hole state, the van Hove singularity is shifted 
towards higher energy part of the band, unoccupied by 
holes. 

Apart from breaking the electron-hole symmetry, the 
extra parameter if modifies the shape of the FS of the free 
electrons and indeed it becomes more consistent with the 
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FIG. 7: Effect of the next-neighbor hopping t'/t = —0.3 
on the noninteracting 2D DOS at the doping x = 1/4: (a) 
hole doping (t = 1); (b) electron doping (t = —1). Dotted 
line shows the Fermi energy in the undoped case, whereas the 
gray area shows the states occupied by either electrons (a) or 
holes (b). 




FIG. 8: FS obtained in the tight binding model at the doping 
x = 1/4: (a) hole doping with t = 1 and: t' = —0.3 (black 
solid line), if = 0.3 (gray solid line), and t' = (dashed line); 
(b) electron doping with t = — 1 and: if = 0.3 (black solid 
line), t' = —0.3 (gray solid line), and t' = (dashed line). The 
long-dashed line in both panels corresponds to the undoped 
case with if — 0. The excessively large value of \t'\ = 0.3 as 
compared to LSCO was chosen only for more clarity of the 
figure. 



FS topology seen by ARPES |M • In the electron- 
doped system NCCO, the low-energy spectral weight at 
the doping x — 0.04 is concentrated in small electron 
pockets around the (±7r, 0) and (0, ±7r) points. Upon 
increasing doping, one observes both the modification of 
the hole pockets and the emergence of new low- lying spec- 
tral weights around (±n/2,±n/2). Finally, at x = 0.15 
the FS pieces evolve into a large holelike curve centered 
at M = (-7T, 7r), In contrast, it has been observed that 
in the lightly doped regime (x = 0.03) doped holes in 
LSCO enter into the hole pockets around (±7r/2, ±7r/2) 
points 0|, implying that the FS is holelike and centered 
at the M point. However, in the heavy overdoped regime 
x = 0.3 it converts into the electronlike FS around the 
T = (0, 0) point. 

Fig. IHa) shows that the model 111211 with if = has a 
nested square FS at half-filling which becomes electron- 
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FIG. 9: Eigenenergy maps of the tight-binding model 11211 
with t' /t = —0.3 as obtained for: (a) hole doping (t = 1); (b) 
electron doping (t = — 1). 

like and shrinks around the T point upon hole doping. 
However, negative t' — —0.3 removes the FS nesting at 
half filling, and the FS expands in the (±fc, 0) and (0, ±k) 
directions, while contracts along the nodal (k, ±fc) and 
(±k,k) directions due to a large gradient dE/dk along 
the latter. Indeed, the eigenenergy map, illustrated in 
Fig. EI a), has in this case a valleylike character with 
a minimum at the T point. Therefore the FS turns 
into a holelike one with experimentally observed arcs [cf. 
Fig. Ufa)]. In contrast, the nearest neighbor hopping t' 
with the same sign as t interchanges the expansion- and 
contraction directions which results in the electronlike 
FS. 

Regarding the electron doped case with t = — 1, shown 
in Fig. IHfb), positive t' = 0.3 (dark solid line) also leads 
to the appearance of arc segments of the FS and makes 
it closer to experimental observations. In this case, how- 
ever, the minimum energy is found at the M point, as 
illustrated in Fig. Efb). It should be noted in passing 
that this FS describes the same situation as the one ob- 
tained with t = 1 and t' — 0.3, indicated by the gray 
solid line in Fig. OJa). In fact, the sign of t is less im- 
portant and turns out to be equivalent to the (ir, ir) shift 
of the momentum without changing the corresponding 
eigenvalues. Consequently, in order to imitate the effect 
of hole and electron doping it is sufficient to study the 
Hamiltonian (QJ only below half-filling and the alterna- 
tion between two regimes is possible by the particle-hole 
transform, 

4 - (-1)*^, (13) 

mapping the model Ill2t with t' < onto the one with 
t' > 0. Therefore, in order to avoid any further confusion 
concerning the signs of t and t' in Eq. (TL2l) , we set here- 
after t to be positive; then a negative t' (t'/t < 0) corre- 
sponds to hole doping, whereas a positive one (t'/t > 0) 
indicates electron doping. 

The remarkable differences of the electronic structure 
due to the broken hole-electron symmetry by t', result 
in different phase diagrams of LSCO and NCCO. In the 
former the long-range AF order is already suppressed in 
the lightly doped regime x ~ 0.03, while in the latter the 
antiferromagnetism is known to be quite robust at in- 



creasing electron doping, hence only commensurate spin 
fluctuations are observed at x — 0.15 [68). The robust- 
ness of the commensurate spin fluctuations in the elec- 
tron doped re gime is consistent with the ED studies of 
the t-t'-J [H |H and t-t'-t"-J E3 models. It is 
also supported by the conclusion that a negative t' pro- 
motes incommensuration at a lower doping level than 
a positive one, reached using the QMC technique ap- 
plied to the extended Hubbard model Finally, the 
XPS measurements in NCCO show that the chemical po- 
tential monotonously increases with electron doping [6J|, 
whereas its shift is suppressed in the underdoped region 
of LSCO jfi6J. These data have been nicely reproduced in 
Ref. jS6J for both compounds, except for the low doping 
regime of LSCO where stripes are expected. All these 
numerical and experimental results indicate that doped 
electrons might selforganize in a different way than holes 
do — in the latter case DWs are formed. Nevertheless, 
stable diagonal stripes with one doped electron per site 
in a DW have been obtained in the slave-boson studies of 
a more realistic extended three-band model [zj, so the 
problem is still open. 

Turning back to the competition between stripes in a 
doped system, Fig. llOf al shows that negative t' stabilizes 
the DSC stripes, whereas positive t' favors the VSC ones, 
within the parameter range where t' does not drive a 
stripe melting. Analogous crossover from vertical stripes 
at small \t'\ to more complex in shape diagonal ones at 
t'/t = -0.1 and t'/t = -0.2 has been found in other HF 
studies The explanation is contained in Ta.hle Fvl 

negative t' gives a positive kinetic energy contribution, 
which is much more readily minimized by the diagonal 
charge configuration. Indeed, despite the solitonic mech- 
anism yielding a noticeable kinetic energy loss due to 
the transverse hopping t'/t = —0.15, the overall kinetic 
energy loss in the case of DSC stripes along the diago- 
nal (11) and antidiagonal (11) directions is smaller than 
the corresponding one for the VSC stripe. A more care- 
ful analysis shows that hole propagation along the DSC 
stripe results in a contribution having the same sign as tl . 
However, it is entirely canceled by the ones coming from 
diagonal bonds of the AF domains so that E*,~ v = 0. 

One observes further that positive t' reduces the 
anisotropy between the kinetic energy gains in the x- and 
y-directions for the VSC stripes, and makes their sum 
more favorable, while negative t' has the opposite effect. 
For the DSC stripes the total kinetic energy also follows 
the same trend. The explanation of these results follows 
from the reinforcement of stripe order by a negative t' (cf. 
values in parenthesis in Table ITTTj) . which suppresses the 
hopping contributions, and its smearing out by positive 
t' where hopping is enhanced. These trends agree with 
the earlier finding within the DMFT that the VSC stripe 
phase is destabilized by kink fluctuations [2jJ. However, 
this stripe (dis) ordering tendency also leads to a consid- 
erably greater change in the Coulomb energy Ejj, listed 
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t'/t 


Ef/t 


E\/t 


E x t r y /t 


E x t + y /t 


Eu/t 


Et t/t 


vsc 


-0.15 


-0.6876 


-0.5886 


0.0140 


0.0140 


0.4778 


-0.7704 


DBC 


-0.15 


-0.6279 


-0.6279 


0.0000 


0.0183 


0.4562 


-0.7813 


DSC 


-0.15 


-0.6275 


-0.6275 


0.0000 


0.0188 


0.4533 


-0.7829 


DBC 


0.15 


-0.6442 


-0.6442 


0.0000 


-0.0282 


0.4883 


-0.8283 


DSC 


0.15 


-0.6437 


-0.6437 


0.0000 


-0.0279 


0.4855 


-0.8298 


VB(S)C 


0.15 


-0.6612 


-0.6372 


-0.0169 


-0.0169 


0.4997 


-0.8325 



TABLE V: Energies per site: ground-state energy Etot, kinetic energy contributions for the bonds along (10) E*, (01) E\ , 
(11) E*~ y and (ll) E®, directions, as well as the potential energy Eu, all normalized per one site, in the extended hopping 
Hubbard model with U/t — 5 and x = 1/8. VBC stripe is unstable at t'/t = —0.15. 




_q 2 I ■ I ■ I ■ I ■ I ■ I -0 2 ' 1 ' 1 ' 1 ' 1 ' 1 

'3.5 4.0 4.5 5.0 5.5 6.0 '3.5 4.0 4.5 5.0 5.5 6.0 

U/t U/t 



FIG. 10: Phase boundaries for: (a) site-centered, and (b) 
bond-centered stripes as obtained in the extended Hubbard 
model with the next-neighbor hopping t' for doping x = 1/8 
(solid line) and x = 1/6 (dashed line). 

in Table El for the DSC than for VSC stripes, which 
contributes significantly to the predominance of the for- 
mer structure for negative t' . In fact, it follows from the 
increase of hole density within the nonmagnetic stripes 
and the magnetization density enhancement within the 
AF domains (cf. Table Hill. 

Like their SC counterparts, DBC stripes are also sta- 
bilized by negative t' resulting in a phase diagram shown 
in Fig. HOf b). In this case, expelling holes from the AF 
domains enhances not only magnetization of their atoms 
but also increases magnetic moment of the hole rich DWs, 
as illustrated in Table HVl This enhancement must, how- 
ever, strongly suppress the dominant transverse kinetic 
energy gain of the VBC stripes. Therefore, the latter are 
already unstable at t'/t — —0.15. 

It is worth noting that a finite diagonal hopping t' 
should directly affect the competition between the d-wave 
pairing correlations and stripes. Indeed, a systematic 
comparison of stripe and pairing instabilities within the 
DMRG framework has shown that when the stripes are 
weakened by positive t' , the latter are strongly enhanced 
due to increasing pair mobility This effect is accom- 
panied by a simultaneous enhancement of the AF corre- 
lations jSjj. Conversely, negative t' reinforcing a static 
stripe order results in the suppression of pair formation in 
the underdoped region, as found both in the DMRG tech- 



nique and Variational Monte Carlo (VMC) 0- How- 
ever, the enhanced pairing correlation, attributed to the 
change of the FS topology in LSCO, has been obtained 
in the optimally doped and overdoped regimes 



Effect of the nearest-neighbor Coulomb interaction 

V 

We now investigate the changes in the stripe sta- 
bility due to either repulsive (V > 0) or attractive 
(V < 0) nearest-neighbor Coulomb interaction, which 
give the phase boundaries between the VSC and DSC 
stripe phases shown in Fig. Illf a). We have found that 
realistic repulsive V favors the latter. The tendency to- 
wards the DSC stripe formation at V > is primarily 
due to a large difference between charge densities at the 
atoms of the DW itself and at all their nearest-neighbor 
sites, a situation which is avoided in the case of VSC 
stripe phases (cf. Fig. EJ. Consequently, the former op- 
timize better the repulsive potential energy component 
Ey, as shown by the data reported in Table PvTl Sim- 
ilarly, the fact that the nearest-neighbor interaction V 
is well minimized only by inhomogeneous charge densi- 
ties makes the DBC stripe phase more favorable than the 
VBC one, as shown in Fig lllf b). While this is also the 




3.5 4.0 4.5 5.0 5.5 6.0 '3.5 4.0 4.5 5.0 5.5 6.0 
U/t U/t 



FIG. 11: Phase diagrams for the site-centered (a) and bond- 
centered (b) stripes obtained in the extended Hubbard model 
with the nearest-neighbor Coulomb interaction V for doping 
x — 1/8 (solid line) and x — 1/6 (dashed line). 
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V/t 


E?/t 


E y Jt 


Eu/t 


E v /t 


Etot/t 


DBC 


-0.4 


-0.6322 


-0.6322 


0.4626 


-0.6194 


-1.4212 


DSC 


-0.4 


-0.6319 


-0.6319 


0.4602 


-0.6193 


-1.4229 


VB(S)C 


-0.4 


-0.6655 


-0.6083 


0.4749 


-0.6251 


-1.4240 


VB(S)C 


0.4 


-0.6838 


-0.6214 


0.5063 


0.6207 


-0.1782 


DBC 


0.4 


-0.6424 


-0.6424 


0.4829 


0.6176 


-0.1843 


DSC 


0.4 


-0.6412 


-0.6412 


0.4789 


0.6171 


-0.1864 



TABLE VI: Energies per site: ground-state energy E tot , ki- 
netic energy (E£,E%) and potential energy (E V ,E V ) com- 
ponents in the extended Hubbard model with the nearest- 
neighbor Coulomb interaction V for U/t — 5 and x = 1/8. 



leading mechanism for both diagonal stripe suppression 
at V < 0, the asymmetry of the curve in Fig. El arises 
from the fact that the lower U values at the transition 
favor the higher kinetic energy contributions available for 
the vertical stripes. 

However, it has been argued based on the results ob- 
tained using the SBA that an increasing repulsive inter- 
action V favors half-filled vertical stripes, hence the latter 
take over at V/t ~ 0.1 in the parameter regime of x = 1/8 
and U/t = 10 (ljj. This finding could naturally explain 
the appearance of filled diagonal stripes in the nickelates, 
provided that they were characterized by a small V term, 
and the stability of the half-filled vertical ones in the Nd- 
codoped cuprates due to possibly larger value of V. It is 
also worth mentioning other HF and variational j^l 
studies in which a variety of intriguing stripe phases, co- 
existing at V/t ~ 1.5 with charge order, has been found 
in a broad doping region. 



Effect of the lattice deformations 

So far, we have demonstrated that a finite anisotropy 
of the transfer integral t can tip the balance between 
vertical and diagonal stripes. Here we will show that 
such anisotropy naturally emerges in a doped system with 
DWs, described by a single-band Peierls-Hubbard Hamil- 
tonian, 

H = - ]T t^u^c^c^ + n iT n H + \ K u l- 

ijv i (ij) 

(14) 

In this model we keep only the leading term and assume 
a linear dependence of the nearest neighbor hopping ele- 
ment on the lattice displacements , 



tij(Uij) =t a (l + au tJ ). 



(15) 



Furthermore, we include the elastic energy oc K which 
allows to investigate the stability of the system with 
respect to a given lattice deformation and to deter- 
mine the equilibrium configuration. For convenience, 
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FIG. 12: Local hole n h (l x ) (top) and magnetization S-n(l x ) 
(second row) density; fractional change of the length for the 
bonds to the right nearest-neighbor along the ^-direction ui°^ 
(circles) and double occupancy D(l x ) (squares) (third row), 
as well as the kinetic energy Ef(l x ) projected on the bonds in 
the z-direction (bottom) of the VSC (left) and DSC (right) 
stripe phases, as obtained in the Peierls-Hubbard model Q14II 
with U/t = 5, A = 0.5 and x = 1/8 (filled symbols). For 
comparison the results obtained with A = are shown by 
open symbols. 



we parametrize the electron-lattice coupling with a sin- 
gle quantity, A = a 2 to/K, with the parameter values 
K/t = 18A" 2 and a = 3A 1 assumed following the 
earlier HF studies [2nj| . As previously, we focus on the 
doping x = 1/8 (a; = 1/6) and present the results of 
calculations performed on 16 x 16 (12 x 12) clusters, re- 
spectively, with periodic boundary conditions. These cal- 
culations have shown that such clusters give the most 
stable filled stripe solutions for the selected doping lev- 
els. The model ifLljl was solved self-consistently in real 
space within the HF (2J. Thereby, we used an approx- 
imate saddle-point formula for the equilibrium relation 
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FIG. 13: The same as in Fig. ^] but for the bond-centered 
stripes. 



between the actual deformation tiy of a given bond and 



the bond-charge density (c\ a c 



,(0) 



^ V(ct c 



■/l.t 



(16) 



being a consequence of the linearity assumption in 
Eq. (HJ. 

Quite generally, it is a widely spread out belief that 
inhomogeneous states at finite doping are very sensitive 
to small changes of A, supported both by the HF 
and ED studies |2f|. Further, it has been shown that 
the electron-lattice interaction favors DW solutions over 
other possible phases, such as isolated polarons or bipo- 
larons [26]. Therefore, a complete discussion of the stripe 
phase stability in correlated oxides has to include the cou- 
pling to the lattice. 

We turn now to the most important aspect of this 
Section. Figs. El an d E3 illustrate the effect of the fi- 
nite electron-lattice coupling A = 0.5 on the SC and BC 



i 


i 


2 


3 4 


5 


VSC (n hi ) 
(St) 


0.270 
0.000 


0.212 
0.146 


0.103 0.038 
0.259 0.310 


0.022 
0.321 


DSC (n w ) 

{s?) 


0.292 
0.000 


0.179 
0.193 


0.094 0.058 
0.277 0.305 


0.046 
0.314 



TABLE VII: Local hole (n hi ) and magnetization ($f ) density 
at nonequivalent atoms of the SC stripe phases, all labeled 
by decreasing hole density in the a;-direction, in the Peierls- 
Hubbard model on a 16 x 16 cluster with U/t = 5, A = 0.5 
and x = 1/8. 



stripes, respectively. Both figures give a clear demonstra- 
tion that, in contrast to the hopping anisotropy ej {JJ 
discussed above, finite A markedly modifies both the lo- 
cal hole density (JHJ and modulated magnetization @ [cf. 
also Table Hill with Ell (SC stripes) and Table HH with 
IVIIII (T3C stripes)]. Basically, the influence of A resem- 
bles the effect of positive t', smearing out the stripe or- 
der by ejecting holes from the DWs, being however much 
stronger. In fact, hole derealization not only suppresses 
the magnetization within the AF domains, but also no- 
ticeably quenches magnetic moments of the BC domain 
walls. These trends can be understood by considering en- 
ergy increments: the kinetic E t , on-site Ey, and elastic 
energy Ek , as explained below. 

One should realize that a system described by the 
Hamiltonian ([14"I) might be unstable towards lattice de- 
formations only if the covalency increase is large enough 
to compensate both the Ejj and Ek energy cost. With- 
out the electron-lattice coupling, a compromise solution 
is mainly reached by developing a strong magnetic or- 
der in the AF domains, where a possible kinetic energy 
gain is irrelevant, and by forming nonmagnetic or weakly 
magnetic DWs with large hole density. As we have al- 
ready shown, transverse charge fluctuations around the 
DWs yield the leading kinetic energy contribution. How- 
ever, enhanced covalency and mixing of the lower ~ td 
and higher ~ + U energy states between a DW and 
the surrounding sites partly delocalize these states and 
increase double occupancy, 



D{l x ) = ("■(z x ,o),T n (**4)),l.)- 



(17) 



i 


i 


2 


3 


4 


VBC (n hi ) 


0.255 
0.074 


0.156 
0.209 


0.063 
0.291 


0.026 
0.319 


DBC (nhi) 

(S!) 


0.248 
0.103 


0.130 
0.243 


0.073 
0.294 


0.049 
0.312 



TABLE VIII: The same as in TableET!|but for the BC stripe 
phases. 
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Indeed, in the A = case, double occupancy D(l x ) 
reaches its maximum at the DWs, as illustrated in 
Figs. El and El The only exception is the DSC stripe 
phase (right panels of Fig. 112)) with the largest D(l x ) in 
the AF domains. As a consequence, the latter is the 
most localized one with the smallest kinetic energy gain 
(cf. Table HJ). 

The situation changes when turning on the electron- 
lattice coupling. When the electrons couple to the lattice 
(A 7^ 0), the bonds contract, and the saddle point values 
of the distortions iflfil) : = (wy) along (10) and (01) 
direction, respectively, are finite. However, a nonuniform 
charge distribution results in a different bondlength in 
the cluster. This is illustrated in Figs . 1 1 21 and 1 1 31 showing 
a fractional change of the length for the bonds to the right 
nearest-neighbor along the ^-direction u x Q ^ (third row). 
Although the values of Uy in the AF domains are also 

substantial, the largest lattice deformations ~ ( c \ a c ja) 
appear either on the bonds connecting atoms of the DWs 
with their nearest neighbors (cf. Fig. 11211 . or on the 
bonds which join two atoms of the bond-centered DWs 
(cf. Fig. Ey . Accordingly, a strengthening nearest neigh- 
bor hopping ltl5|l enables a larger kinetic energy gain on 
these bonds (cf. bottom of Figs . IT2l and 113)) . 

As expected, the increasing covalency is accompanied 
by partial quenching of magnetic moments. In order to 
appreciate this tendency, let us consider a site in the 
AF domain with larger density of t-spin electrons (at 
A sublattice). Once the magnetization is reduced, the 
corresponding T-spin energy level which belongs to the 
lower Hubbard band is pushed upwards, and the t-spin 
of the upper Hubbard band goes down. As a result, the 
locally raised t-spin state becomes stronger mixed with 
J,-spin states at the surrounding sites of B sublattice, and 
simultaneously bond-charge density increases. At the 
same time, electrons, jumping forth and back between 
the central site with the t-spin polarization and its near- 
est neighbors with the t-spin one, enhance considerably 
double occupancy D(l x ), as shown in Figs. \12\ and 1131 
This weakens the stripe order and results in a more uni- 
form distribution of D(l x ). Of course, the increase of the 
elastic energy and concomitant enhancement of the on- 
site energy, both owing to finite bond contractions (IKill . 





Ef/t 


EZ/t 


Eu/t 


E K /t 


Etot/t 


DBC 


-0.9679 


-0.9679 


0.6478 


0.2548 


-1.0332 


DSC 


-0.9670 


-0.9670 


0.6450 


0.2544 


-1.0346 


VB(S)C 


-1.0496 


-0.9248 


0.6719 


0.2638 


-1.0387 



TABLE IX: Ground-state energy E to t per site, kinetic energy 
(Ef,Ef) and potential energy (Eu,Ek) components, as ob- 
tained in the Peierls-Hubbard model. Parameters: U/t = 5, 
A = 0.5, and x = 1/8. 
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U/t U/t 

FIG. 14: Phase diagrams for site-centered (a) and bond- 
centered (b) stripe structures as calculated from the Peierls- 
Hubbard model for doping x = 1/8 (solid line) and x = 1/6 
(dashed line). 



is compensated by the kinetic energy gain and the total 
energy is lowered (cf. Tables III and IIXII . 

We close this Section with the phase diagrams shown 
in Fig. El They were obtained by varying U and the 
coefficient a, while maintaining constant K/to = 18A -2 . 
The increased stability of vertical stripes follows from the 
relative stronger enhancement of the local hopping ele- 
ments II 1511 (and consequently larger gain of the kinetic 
energy), especially on the bonds in the direction perpen- 
dicular to the DWs itself. 



SUMMARY 

In summary, we have shown that a competition be- 
tween magnetic energy of interacting almost localized 
electrons and the kinetic energy of holes created by 
doping leads to the formation of new type of coexist- 
ing charge and spin order — the stripe phases. We 
have shown that vertical (horizontal) and diagonal stripes 
dominate the behavior of the charge structures formed by 
doping the antiferromagnet away from half filling, using 
the solutions obtained for the Hubbard model within the 
HF approximation in the physically interesting regime 
of the Coulomb interaction. The detailed charge distri- 
bution and the type of stripe order depend on the ratio 
U/t, on the value of the next-neighbor hopping t', and 
on the nearest-neighbor Coulomb interaction V. We have 
also shown that a strong electron-lattice coupling might 
be responsible for the appearance of the vertical stripes 
observed in the superconducting cuprates at x = 1/8. 

Altogether, although some experimentally observed 
trends could be reproduced already in the HF approach, 
the presented results indicate that strong electron corre- 
lations play a crucial role in the stripe phases and have 
to be included for a more quantitative analysis. Fur- 
ther progress both in the experiment and in the theory is 
necessary to establish the possible role of stripes in the 
phenomenon of high temperature superconductivity. 
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